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Elastic and inelastic neutron scattering studies of Li(Nii_ :E Fe a: )P04 single crystals reveal anoma- 
lous spin-wave dispersions along the crystallographic direction parallel to the characteristic wave 
vector of the magnetic incommensurate phase. The anomalous spin- wave dispersion (magnetic soft 
mode) indicates the instability of the Ising-like ground state that eventually evolves into the incom- 
mensurate phase as the temperature is raised. The pure LiNiPCU system (x — 0), undergoes a 
first-order magnetic phase transition from a long-range incommensurate phase to an antiferromag- 
netic ground state at Tn = 20.8 K. At 20% Fe concentrations, although the AFM ground state is 
to a large extent preserved as that of the pure system, the phase transition is second-order, and the 
incommensurate phase is completely suppressed. Analysis of the dispersion curves using a Heisen- 
berg spin Hamiltonian that includes inter- and in-plane nearest and next-nearest neighbor couplings 
reveals frustration due to strong competing interactions between nearest- and a next-nearest neigh- 
bor site, consistent with the observed incommensurate structure. The Fe substitution only slightly 
lowers the extent of the frustration, sufficient to suppress the IC phase. An energy gap in the dis- 
persion curves gradually decreases with the increase of Fe content from ~2 meV for the pure system 
= 0) to ~0.9 meV for x = 0.2. 

PACS numbers: 75.25.+Z, 75.50.Ee, 78.20.Ls 



I. INTRODUCTION 

Spontaneously occurring incommensurate (IC) struc- 
tures can be classified into two general categories. The 
first group consists of systems for which the IC phase is 
the ground state, and the second group encompasses sys- 
tems for which the IC phase manifests itself as an inter- 
mediate state between a commensurate ground state and 
a highly symmetric phase at higher temperatures 
Systems with incompatible interactions among nearest 
and next nearest neighbors that may lead to geomet- 
rical frustration, in general belong to the first group 
settling into an IC ground state Similarly, nearest- 
neighbors frustrations brought about by off-diagonal 
Dzyaloshinskii-Moriya type interactions, that compete 
with the isotropic interactions, can also give rise to IC 
ground states Q. Magnetic systems consisting of in- 
teracting localized moments, such as MnSi[f|, FeGeJ(|, 
NiBr 2 0, Ba 2 Cu 2 Ge 2 7 [l], CuB 2 4 ,i, LiCuV0 4 [l3, 
and CdCr 2 04[llI are typical examples of the first group. 
On the other hand, the intermediate IC phases are in 
general electronically driven by instabilities due to the 
incompatibility in the interactions of a collective mode 
(phonon) and conduction electrons at the Fermi sur- 
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face. These structurally modulated phases are gen- 
erally observed in metallic systems as charge density 
waves [l2l[l3j. or martensitic transitions in alloys 0, [l^, 
and occur at intermediate temperatures between a disor- 
dered state at high temperatures and a highly symmetric 
ground state at low temperatures |3J. Systems belong- 
ing to this second group possess a few distinct charac- 
teristics: 1) they undergo a first-order commensurate- 
incommensurate (C-IC) phase transition, 2) they give 
rise to strong diffuse scattering above and below the C-IC 
transition 3) they exhibit anomalies in their phonon dis- 
persion curves that signal the emergence of the IC phase. 
A typical phonon anomaly appears as a minimum, or a 
dip, in the dispersion curve, commonly referred to as a 
soft-mode, at a wave-vector that defines the propagation 
vector and the shortest wavelength of the IC modula- 
tion. Due to the first-order nature of the C-IC transi- 
tion, phonons are not well defined close to the transition, 
thus the whole dispersion curve, including the soft mode, 
abruptly disappears near the transition, and in turn, a 
frozen phonon sets in giving rise to a single peak at en- 
ergy uj « 0. The frozen phonon, realized as an elastic or 
quasi-elastic peak at and around the wavevector defining 
the IC structure, is identified with the IC structure. 

An intriguing IC magnetic phase, with features that 
characterize the second group, has been found recently 
in the magnetoelectric crystal LiNiP04. The IC phase 
occurs over a narrow range of intermediate temperatures 
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between an antiferromagnetic (AFM) ground state and 
a high temperature paramagnetic phase [I6fl. Here, it 
was found that LiNiP04 undergoes a first order tran- 
sition from the antiferromagnetic ground state to a long- 
range IC order at a Neel temperature, T/v = 20.8 K 
(T/v = Tc-ic)- As the temperature is increased, a 
second-order phase transition from long-range incom- 
mensurate magnetic order to the paramagnetic state oc- 
curs at Tic = 21.7 KpH fl7l |. The incommensurate 
spin correlations gradually weaken and the spins are es- 
sentially uncorrelated by T rs 35 K. In addition to ex- 
hibiting a first order C-IC phase transition, strong dif- 
fuse scattering below and above the transition has also 
been reported [lg|. This unusual magnetic intermediate 
IC phase has characteristics that classify it with the sec- 
ond group mentioned above, however, it should be noted 
that LiNiPC>4 is an insulator (with an energy gap of ap- 
proximately 1 eV), thus the IC phase cannot be induced 
by interaction with conduction electrons. 

A recent neutron scattering study [18j investigated the 
spin dynamics of pure LiNiP04 to determine the spin 
Hamiltonian and identify other features that character- 
ize the aforementioned second group of IC systems, par- 
ticularly looking for a soft-magnetic mode, the analog of 
the soft mode in structurally IC systems. An unusual 
minimum in the spin-wave dispersions in the AFM com- 
mensurate ground state was observed at the modulation 
vector of the IC phase, and was explained as the precur- 
sor of the C-IC phase transition that originates from a 
trade off between competing Heisenberg interactions of 
nearest and next-nearest neighboring Ni 2+ ions and an 
extra lock-in energy at lower temperatures originating 
from the strong single ion anisotropics found in the sys- 
tem pl| . We have recently reported on the spin dynam- 
ics and magnetic properties of the isostructural LiFeP04, 
LiCoP04 and LiMnP04 systems and found no evidence 
for an IC phase and no anomalous spin- wave dispersions 
[T^ . l2(il . l21j ]. In the present study we have substituted 
Fe for Ni to form LiNii_ x Fe :E P04 single crystals with 
up to x — 0.2 to compare with the magnetic behavior 
of the parent material and shed further light on the IC 
phase in LiNiP04. Our studies show that, up to a sub- 
stitution level greater than ~ 0.15 iron, the IC phase is 
still present, and only at higher Fe concentrations does 
it disappear completelvp2l|. We report herein the spin 
dynamics of LiNio.gFeo.2PO4, that do not show evidence 
for the IC magnetic structure, and compare the results 
with measurements of pure LiNiP04. 

LiNiP04 is an insulator belonging to the olivine fam- 
ily of lithium orthophosphates L1MPO4 (M — Mn, Fe, 
Co, and Ni) with space group Pnma[23j. All mem- 
bers of this family were found to be antiferromagnets 
with the same magnetic structure differing only in the 
spin-direction [24], [2f| Hfjj], however a recent single crystal 
study of LiNiP04 revealed that the magnetic spins are 
not co-linear in the AFM ground state but are slightly 
canted within the ac-plane [27| . Neutron scattering 
studies demonstrated that LiMP0 4 (M =Ni,Co,Mn) ex- 



hibit properties between two-dimensional (2D) and three- 
dimensional (3D) with an interlayer coupling that is 
stronger relative to the coupling found in the cuprates, 
for instance fl9L I20I l2ll l28l 123] . These insulators also ex- 
hibit strong linear magnetoelectric (ME) effects, with the 
observed ME tensor components, a xz , a zx , for LiNiP04, 
in agreement with the antiferromagnetic point groups 
mm'm, but with some anomalies [3Ct |3 lL |32| |. In par- 
ticular, the ME effect measurements of LiNiP04 as a 
function of temperature reveal a first-order AFM tran- 
sition, and an unusual decrease of the ME coefficient at 
temperatures below a maximum close to Tv (l6l |30T| ■ Re- 
cently, a microscopic model combining super- exchange 
and Dzyaloshinsky-Moriya interactions with elastic dis- 
placements of exchange mediating ions has accurately 
explained the temperature dependence of the ME coeffi- 
cients in LiNiP04 (27[. The model shows that the sharp 
decrease of a xz and a zx as a function of temperature is 
intimately connected to the first-order nature of the C-IC 
phase transition in LiNiP04. By contrast, the isostruc- 
tural materials, LiCoP04, LiFeP04, and LiMnPO_4, all 
exhibit a continuous change of the ME coefficients [3l[ re- 
flecting the second-order nature of their magnetic phase 
transition from a commensurate AFM state to the para- 
magnetic state. Magnetic susceptibility studies of poly- 
crystallinc LiNiP04 showed a significant deviation from 
the Curie- Weiss law in a temperature range much higher 
than T/v, and neutron scattering from the same poly- 
crystalline sample gave rise to diffuse scattering at the 
nominal position of the AFM Bragg reflection up to 
T w 2TW[28l|. Recent magnetic susceptibility measure- 
ments of single crystal LiNiP04 showed two features, one 
at T^ = 20.8 K and one at Tic — 21.7 K associated 
with an AFM transition and an intermediate IC phase 
[l7| . in agreement with the observed neutron diffraction 
data p]|. 



II. EXPERIMENTAL DETAILS 

LiNiP04 single crystals were grown by the stan- 
dard flux growth method (LiCl was used as the 
flux) from a stoichiometric mixture of high purity 
NiCl 2 (99.999%) and Li 3 P0 4 (99.999%) H. To pre- 
pare LiNio.8Feo.2PO4, the Fe substitution was intro- 
duced by adding FeCl 2 (99.999%) to the flux at a 
molar ratio of 1:4 to MCI2. The composition of 
Li(Nio.8Feo.2)P04 single crystals were confirmed by 
chemical analysis. X-ray diffraction and GSAS refine- 
ment show that Li(Nio.8Feo. 2 )P04 has the same crystal 
structure and symmetry group as pure LiNiP04. 

The magnetic susceptibility measurements were per- 
formed on a superconducting quantum interference de- 
vice (SQUID) magnetometer. The single crystals used 
for magnetic property measurements were oriented using 
Laue back scattering x-ray diffraction. For the different 
measurements, the single crystals were glued to a plas- 
tic straw with the specified axis parallel to the applied 
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magnetic fields with error less than 5°. 

Elastic and inelastic neutron scattering studies of 
LiNiPC>4 were performed on the HB1A spectrometer at 
High Flux Isotope Reactor (HFIR) at Oak Ridge Na- 
tional Laboratory. A monochromatic neutron beam of 
wavelength A = 2.37 A (14.61 meV, k a = 2tt/\ = 
2.66A- 1 ) was selected by a double monochromator sys- 
tem, using the (002) Bragg reflection of highly oriented 
( mosaicity 0.3 deg) pyrolytic graphite (HOPG) crys- 
tals. The collimating configuration 40'-40'-Sample-34'- 
68' was used throughout the experiments, yielding an 
average energy resolution of ~ 1 meV. Two sets of HOPG 
crystals, located between and after the monochromator 
crystals, were used as filters removing the A/2 compo- 
nent from the incident beam to better than one part in 
10 4 . Elastic neutron scattering from Li(Ni .sFe .2)PO4 
single crystals was measured on the HB1A spectrometer, 
and the inelastic neutron scattering was measured on the 
Spin Polarized Inelastic Neutron Spectrometer (SPINS) 
at the NIST Center for Neutron Research (NCNR) using 
a fixed final energy of 5 meV. The collimating configura- 
tion 80'-Sample-Be filter-80' was used for these measure- 
ments yielding an energy resolution ss 0.2 meV. 



III. EXPERIMENTAL RESULTS 

A. Magnetic Susceptibility 

Magnetic susceptibility measurements of LiNiP04 and 
Li(Ni .8Feo.2)P04 single crystals along the easy c-axis 
are shown in Fig. [1] (a). The measured magnetic suscep- 
tibilities of the two systems are different in two respects. 
First, the absolute value of the susceptibility in Fig.[T](a) 
is larger for the iron substituted sample, indicating the 
presence of uncompensated paramagnetic sites, due to 
the random distribution of iron spins (S = 2) with a mo- 
ment that is different than that of Ni +2 (S = 1). We note 
that susceptibility measurements under field- or zero-field 
cooling indicate subtle spin-glass properties [Hj]. Second, 
large differences are identified in the derivatives of the 
susceptibilities with respect to temperature, as shown in 
Fig. [T] (b ). I n agreement with previous measurements of 
LiNiP04[17|, the main AFM-IC transition has the char- 
acteristics of a first order transition, and the anomaly as- 
sociated with the transition from the long-range IC struc- 
ture to the the paramagnetic state at Tjc = 21.7 K is in 
good agreement with the neutron diffraction studies [l6l|. 
By contrast the Fe substituted (x = 0.2) crystal has only 
one smooth feature characteristic of a second order phase 
transition with no indication of a secondary transition. 
We therefore conclude that whereas the Fe substitution 
maintains the AFM ground state for all x < 0.2 it does 
not modify the nature of the transition and does not elim- 
inate the IC long-range order up to a substitution level 
of z~ 0.2. 
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FIG. 1: (color online) (a) Susceptibility measurements of 
LiNiP04 (solid line) and Li(Ni .8Fe .2)NiPO4 (dashed line), 
(b) The respective derivatives of the susceptibilities with re- 
spect to temperature, showing the two features due to the 
transitions to IC and AFM in pure LiNiP04 and a single 
broad feature (second order transition) in Li(Nio.sFeo.2)P04 



B. Elastic Neutron Scattering 

Neutron diffraction measurements confirm the or- 
thorhombic structure of both the pure and the Fe sub- 
stituted samples. For pure LiNiP04 we find the follow- 
ing room temperature lattice parameters of a — 10.030, 
b = 5.847, and c = 4.677 A and for Li(Ni . 8 Fe .2)PO 4 
we find a = 10.057, b = 5.881, and c = 4.672 A at 10 
K. For the x — 0.2 sample, we identified a weak nuclear 
peak (at room temperature), not identified in the X-ray 
diffraction of the powder, at the (010) that may indicate 
a small structural distortion along the 6— axis. In gen- 
eral, elastic neutron scattering and magnetic susceptibil- 
ity measurements of the Fe substituted Li(Nii_ a; Fe :i ;)P04 
crystals show the low-temperature ground states of these 
systems are antiferromagnetic with a magnetic arrange- 
ment similar to that found in pure LiNiP04. 

The magnetic spins in LiNiP04 are primarily directed 
along the c-axis in the AFM ground state, but are slightly 
canted with a small component along the a-axis. As the 
Fe concentration is increased, up to at least x w 0.2, 
the Neel temperature changes slightly but the nature of 
the order parameter changes more dramatically. Figure 
[5] shows the temperature dependencies of the magnetic 
order parameters for LiNiP04 and Li(Nio.sFeo.2)P04 as 
measured on the (010) magnetic peak. As previously 
discussed [l6|. LiNiP04 undergoes a first order magnetic 
phase transition from commensurate AFM ground state 
(labeled A in Fig. ^ to a long-range incommensurate 
structure at T/v = 20.8 K, and subsequently to the para- 
magnetic state at Tic ~ 21.7 K by a second-order phase 
transition. The incommensurate spin correlations are 
gradually lost by a temperature between 34 K to 40 K. In 
contrast, the transition from the AFM to the paramag- 
netic phase in Li(Nio.sFeo.2)P04 is continuous, i.e., it is a 
second-order phase transition to the paramagnetic phase 
with no clear evidence for any intermediate magnetic 
phases. The temperature dependent order-parameter for 
Li(Nio.8Feo.2)P04, shown in Fig. [21 was fit to a power- 
law function (solid line) yielding a transition temperature 
T N = 20.6±0.2 K and a critical exponent = 0.33±0.03. 
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FIG. 2: (color online) Magnetic order parameter of 
pure LiNiPCU and Li(Nio.sFeo.2)P04 versus temperature. 
LiNiPCU undergoes a first order phase-transition from antifer- 
romagnetic ground state to long-range incommensurate struc- 
ture at Tiv = 20.8 K, and at T IC » 21.7 K the IC structure 
transforms to the paramagnetic state. The incommensurate 
spin correlations become negligible at about 35 to 40 K. By 
comparison, Li(Nio.8Feo.2)P04 transforms from the collinear 
ground state to the paramagnetic state by a second order 
phase transition at Tjv = 20.6 K. The labeled temperature 
regions refer to phases of pure LiNiPCU. 

Figure [3] shows scans along the (OkO) direction for the 
pure and x = 0.2 samples above the Neel temperature. 
Whereas these scans for LiNiP04 above 20.8 K show two 
satellite peaks due to the IC phase, with intensities, peak- 
shapes, and wave-vectors that are strongly temperature 
dependent, no similar peaks along the (OkO) direction 
or along any other principal direction were detected for 
Li(Nio. 8 Feo.2)P0 4 . 



C. Inelastic Neutron Scattering 

Figure HJa) shows constant-Q energy scans of spin- 
waves propagating along (OqO) for LiNiP04 at 10 K mea- 
sured on the HB1A spectrometer at HFIR (energy resolu- 
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FIG. 3: (color online) Longitudinal scans along the (OkO) di- 
rection for (a) Li(Ni .8Fe .2)PO4 and (b) L1MPO4. Above the 
Neel temperature no satellite peaks, due to the long-range in- 
commensurate structure, were observed for LiNio.sFeo.2PO4 
in the temperature range between 18 K to 22 K with 0.25 K 
temperature steps. The error bars in this paper are statistical 
in origin and represent one standard deviation, (r.l.u, stands 
for reciprocal lattice units, for example for the (OqO) direction 
q is normalized to b* — 2ir/b.) 



tion ss 1 meV). Similar constant-Q energy scans obtained 
on the SPINS spectrometer at NCNR (energy resolution 
0.2 meV at zero energy transfer) on Li(Nio.gFeo.2)P0 4 
at 4 K are shown in Fig. Efb). Each constant-Q scan 
was fit to a Gaussian profile (including a constant back- 
ground) shown as solid lines. Using this analysis, the 
spin-wave dispersion curves along all the three princi- 
pal directions, (£00) (0£0) and (00£) for both crystals 
were compiled in Fig. [5] The dispersion curves show 
an energy gap that decreases with iron substitution. A 
gap of AE ~ 1.9 meV is observed for LiNiPO^ com- 
pared with AE ~ 0.9 meV for Li(Nio. 8 Fe .2)P0 4 . The 
dispersion curves along the propagation vector (0g0) of 
the AFM structure are softer (lower in energy) than the 
curves along the other principal directions. In particu- 
lar, it is even softer than inter-layer spin-waves, along 
the (qOO) direction, propagating perpendicular to the b- 
c planes. This behavior should be contrasted with the 
spin-waves of isostructural LiFcP04, where the disper- 
sion along the (OqO) direction (l9j is stiffer than that along 
the (qOO). Most importantly, for small g's the curve is 
almost flat, with a shallow minimum at q w 0.1, i.e., a 
soft magnetic mode, whereas for a simple gapless AFM 
systems the spin-wave dispersion is expected to be lin- 
ear at small wave-vectors. We identify the anomalous 
spin-wave dispersion along (0,g,0) direction with the soft 
magnetic mode. 
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FIG. 4: (color online) (a) Constant-Q energy scans along 
the (OqO) for a single crystal LiNiP04 at 10 K (b) and for 
Li(Nio.8Feo.2)P04 at 4 K. The solid lines are Guassian fits in- 
cluding constant background. All modes were measured from 
the (010) zone center. 



The substitution of Fe in Li(Ni 1 _ a; Fe a ;)P0 4 modifies 
both the energy gap and the overall dispersion curves. In 
particular, the mode along the (OqO) direction is modified 
and the shallow minimum is not observed, as shown in 
Fig. H]for Li(Nio.8Feo.2)P0 4 . These modifications in the 
spin- wave dispersion, are not striking considering the fact 
that the IC phase is not present in this crystal. This may 
suggest that although the ingredients for the IC phase 
to occur are still present, namely competing interactions 
that lead to frustration, they are not sufficiently strong 
or coherent to stabilize an equilibrium IC phase above 
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FIG. 5: (color online) Spin wave dispersion curves along the 
(OOq), (0g0) and {q00) directions at (a) 10 K for LiNiP0 4 , (b) 
at 4 K for Li(Nio.8Feo.2)PC>4. The solid and dashed lines are 
fits using the spin- wave Eq. @. The dotted line starting at 
the zone center indicates the spin wave optical branch. 



IV. ANALYSIS AND DISCUSSION 




FIG. 6: (color online) (a) Atomic structure of LiNiP04. The 
magnetic moments of Ni 2+ are along c-axis in L1MPO4. (b) 
Illustration of spin couplings in LiNiP04. The same defini- 
tions of bonding are used in the spin wave Hamiltonian 



A. Linear Spin- Wave Theory 

To anal yze the measured spin-wave data we follow the 
model of 1811 using linear Holstein-Primakoff spin- wave 
theory [lg, [34| to calculate the eigenvalues as a function 
of wavevectors. The interaction parameters determining 
the eigenvalues are then refined by a non-linear-least- 
square fit to the measured dispersion curves. LiNiP04 
adopts the Puma symmetry group, in which Ni 2+ (5 = 1) 
ions occupy the centers of slightly distorted NiC>6 oc- 
tahedra, and P ions are located at the centers of PO4 
tetrahedra. The NiC>6 octahedra are corner shared and 
cross-linked with the PO4 tetrahedra forming a buckled 
two-dimensional plane normal to the a-axis. The atomic 
structure and definition of spin coupling, used in this 
study, are illustrated in Fig. [6j The small canting of the 
Ni 2+ s pins has a negligible influence on the spin wave 
model 18] so for simplicity we have assumed a ground 
state with spins pointing strictly along the c axis. The in- 
plane nearest-neighbor (NN) coupling ( Ji) is mediated by 
an oxygen through a Ni 2+ -0-Ni 2+ bond. The distances 
between the inplane NN are 3.806 A. There are two in- 
plane next-nearest-neighbors (NNN), with distances of 
5.89lA and 4.705A, with inplane couplings J 2 and J3, 
respectively. These NNN are linked via Ni 2+ -0-P-0- 
Ni 2+ bond. For inter-layer coupling, we consider only 
the NN interactions J4 and J5 in adjacent layers (5.397 
and 5.495A apart, respectively). The exchange interac- 
tion between NN in adjacent layers is through phosphate 
tetrahedra. The spin-coupling via phosphate tetrahedra 
can be significant and cannot be ignored, as has been 
found for Li 3 Fe 2 (P04)3_where all spins are coupled via 
phosphate tetrahedra [35j. In addition to the Heisen- 
berg interactions, the spin Hamiltonian includes standard 
single-ion anisotropy terms D^(S^) 2 (£ l = x,y,z,) as fol- 
lows, 



H = ^(J {!j} S l .S J ) + ^^(5«) 2 



(1) 



where D XjVtZ are the single ion anisotropics along the a, 
b and c axis respectively. Since the excitation spectrum 
is insensitive towards an overall shift in the ground state 
energy we define D z = for simplicity. The magnon dis- 
persion derived from Eq. |T]) by linear spin-wave theory 
is given by Eq. ([2]). 



hw = \J A 2 — (B ± C) 2 , 



(2) 



where, 
A = 



45(Ji + J B ) - 25[J 2 (1 - cos(q • r 5 )) + J 3 (l - 
cos(q • r 6 )) + J 4 (2 - cos(q • r 7 ) - cos(q • r 8 ))] 
+ (S-l/2)(D x + D y ), (3) 



B = (S-l/2)(D x -D y ), 



(4) 



C= 2J 1 5[cos(q • ri) + cos(q • r 2 )] 

+2J 5 5[cos(q • r 3 ) + cos(q • r 4 )], (5) 

and Yi denotes a vector to a NN and NNN, ri = 
(0,6/2, c/2); r 2 = (0, 6/2, -c/2); r 3 - (a/2, 6/2, 0); r 4 = 
(a/2,-6/2,0);r 5 = (0,6,0);r 6 = (0,0,c);r 7 = 
(a/2, 0, c/2); r 8 = (a/2, 0, —c/2). In our model, the calcu- 
lated spin waves have two non-degenerate branches (de- 
noted by the ± sign in Eq. \2§ as a result of the different 
anisotropies along the x, y and z direction. 

The energy gaps at q = for the two branches are 



AE = JWS(S - l/2)U B (Ji + J 5 ) + 4(5 - l/2) 2 D x D y , 



for (B - C) in equation [2] and 



(6) 



AE = ^16S(S - l/2)D y (J 1 + J 5 ) + 4(5 - l/2) 2 D x D y , 

(7) 

for (B + C). The equations show that the energy gaps 
depending on both the single-ion anisotropy terms and 
the exchange interactions. 
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5 = 1 and S — 1.2 are used for LiNiPC>4 and 
Li(Ni .sFe .2)PO4, respectively. The experimental data 
for L1N1PO4 and for Li(Nio.sFeo.2)P04 were simultane- 
ously fit for the three principal directions by Eq. [5] using 
the (B - C) dispersion. The best fits, shown by solid 
lines in Fig. [5] were obtained by using the parameters 
listed in Table Q] It is noted that the values given in 
Table U for pure LiNiPC>4 are consistent with those re- 
ported in Ref. [l8|. The dashed lines in Fig. [5] are the 
second mode of the spin wave calculated using (B + C) 
in Eq. [5] and the parameters listed in Table fl] For the 
spin wave dispersion along the (0, 1 + q, 0) direction, sev- 
eral excitations for the second mode were measured. It 
is clearly shown in Table Q] that the inplane NN exchange 
Ji is much larger than the inter-plane NN exchanges, J4 
and J5, consistent with the the quasi-2D character of the 
system. The coupling constants also show that the NNN 
inplane coupling along the 6-axis J2 has the same sign 
as that of J 1; implying competing interactions. In par- 
ticular wc find that J2, which couples spins along the 
6-axis, is significantly larger than J 3 that couples NNN 
along the c-axis. This is the direction along which the IC 
structure is realized. The single ion anisotropies, D x and 
D y , are both positive indicating that a c-axis magnetic 
moment is a favorable ground state, as observed experi- 
mentally. The Fe substitution systematically weakens all 
effective spin-couplings and the single ion anisotropies. 
Two non-degenerate branches of the spin wave disper- 
sion have been observed at several scattering vectors in 
LiNiP04 using the high-flux thermal neutron triple axis 
IN8 at Institut Laue-Langevin (ILL) and were perfectly 
fit by the proposed spin- wave model [l8j . 



TABLE I: Best fit exchange parameters and single ion 
anisotropies used to fit the dispersion curves in Fig. [5] 



LiNiP0 4 



LiNio.gFeo.2PO4 



Ji 
Ji 
Js 
Ji 
Js 
D x 

Dy 

D z 



0.94(08) 

0.59(05) 

-0.11(05) 

-0.16(02) 

0.26(02) 

0.34(06) 

1.92(01) 





0.88(15) 

0.44(04) 

0.087(02) 

-0.22(04) 

0.038(004) 

0.072(006) 

1.47(1) 





V. SUMMARY 

Model calculations of spin systems with competing in- 
teractions between NN and NNN have demonstrated that 
anomalous spin- waves, i.e., soft-maqnetic-mode are pos- 
sible for such frustrated systems [3(1 [37| • The spin cou- 
plings for the Fe substituted compound [x — 0.2) are 
slightly different than those of the the pure one with sim- 



ilar frustrations, but they do not lead to the IC phase. 
The realization of the IC phase as an intermediate state 
may be related to the energy gap compared to thermal 
energies at Tjv. It is interesting to note that the energy 
gap observed in the dispersion curves of the pure system 
is very close to k B T N , i.e, AE w 1.9 meV = 22 K. By 
contrast the energy gap in the Fe substituted system is 
much lower than the intrinsic T/v temperature AE w 0.9 
meV = 10 K. Thus, although the ingredients for the IC 
phase are present in the Fe substituted sample and give 
rise to diffuse scattering, they cannot stabilize the IC 
structure at any temperature. Another measure for the 
feasibility of an IC phase is the ratio of the competing in- 
plane couplings J2/J1, for example. The reduction of the 
ratio from J 2 /Ji « 0.63 for LiNiP0 4 to J2/J1 « 0.5 for 
the Fe substituted system (x = 0.2) is sufficient to desta- 
bilize the IC phase. LiFeP0 4 LiFeP04, with J 2 /Ji ~ 0.4, 
exhibits a second-order paramagnetic-AFM phase tran- 
sition with no evidence for the IC magnetic structure at 
any temperature . 

In summary inelastic neutron scattering studies of 
LiNiP04 and Li(Nio.sFeo.2)P04 show the spin-dynamics 
of these systems is anomalous. Whereas the anomaly in 
the pure material leads to an IC intermediate state, the 
reduced anomaly in the perturbed system with the sub- 
stitution of Fe for Ni does not exhibit an IC magnetic 
structure. The spin wave dispersion curves for both sys- 
tems were analyzed using the eigenvalues obtained from 
a Heisenbcrg-like spin Hamiltonian by linear spin wave 
theory. The spin couplings obtained indicate frustration 
between inplane NN and NNN, in particular along the 
direction that the IC structure is observed. Although Fe 
substitution does alter the ground state, and preserves 
the frustration to a lesser degree, it eliminates the IC 
phase altogether. 
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